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We have constructed a GEANT4-based detailed software model of photon transport 
in plastic scintillator blocks and have used it to study the NEMO-3 and SuperNEMO 
calorimeters employed in experiments designed to search for neutrinoless double beta 
decay. We compare our simulations to measurements using conversion electrons from 
a calibration source of 207 Bi and show that the agreement is improved if wavelength- 
dependent properties of the calorimeter are taken into account. In this article, we briefly 
describe our modeling approach and results of our studies. 

1. Introduction 

Several dedicated efforts have recently been proposed to describe optical photon trans- 
port in scintillator detectors using various Monte Carlo packages [ [lj [2j [3l HJ E] . Using the 
GEANT4 (version 4.9.1 patch 3) framework [[6], we have constructed a comprehensive 
and detailed model of photon transport in plastic scintillator, then used the model to 
study the individual NEMO-3 calorimeter modules. In the model, we account for the 
wavelength dependence of optical properties of the scintillators, light guides, reflective 
wrappings, photodetectors and coupling materials. We use wavelength dependent self- 
absorption and re-emission in the scintillator and account for the fluorescent quantum 
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yield of the wavelength shifter. Our results show that this detailed modeling exhibits 
better agreement with measurements compared to a monochromatic approach. 

The NEMO-3 experiment, located at the Laboratoire Souterrain de Modane in the 
Frejus tunnel, searches for neutrinoless double beta decay by employing tracking and 
calorimetry systems and has been taking data since 2003 [0 EJ [9j [TOj |TTJ. The calorimeter 
modules consist of large polystyrene scintillator blocks with light guides coupled to either 
flat or hemispherical photomultiplier tubes (PMTs). Signals in an individual block are 
due to incident particles, mostly /3 and 7 rays, and the response varies with the energy 
and the impact point on the entrance face. The response also depends on the size and 
geometry of the blocks. The energy resolution and background rejection improves if a 
correction for the non-uniformity due to the impact position is applied for each type 
of employed blocks [[?]. We have reproduced the spatial dependence of response of the 
NEMO-3 scintillators to 207 Bi conversion electrons and have optimized the new scintillator 
block geometry for the next generation double beta decay experiment, SuperNEMO. 



2. Modeling details 
2.1. The detector 

NEMO-3 calorimeter modules [ [7] were manufactured to conform to the overall cylin- 
drical geometry of the detector. Each module faces the isotopic foil (a source of double 
beta transitions) and is composed of a scintillator block, a light guide, and a 3 inch or 
5 inch PMT. The detector is azimuthally divided into 20 identical wedge sectors, each 
assembled as a tracker and calorimeter with a source foil, as depicted in Figure [TJ The 
scintillators hermetically cover the two cylindrical walls which surround the foil and the 
tracking volume. There are three types of blocks comprising the walls and there is also 
partial scintillator coverage of the top and bottom end-caps. On the interior wall, there 
are two IN blocks which are mirror-symmetric. On the exterior wall, there are two types 
of blocks: center (EC) and edge (EE) type blocks, with EE blocks symmetrically placed 
on either side of an EC block. The dimensions of these scintillator blocks are given in Ta- 
ble [TJ Simulations of the response of these three types of blocks was the primary objective 
of our work reported here. 

The scintillator material is composed of polystyrene (PS) (98.49% by weight), a pri- 
mary dopant p-terphenyl (pTP) (1.5% by weight) and a wavelength shifter l,4-bis(5- 
phenyloxazol-2-yl) benzene (POPOP) (0.01% by weight). The light guides are made out 
of polymethylmethacrylate (PMMA) and serve as the interface between the PMTs and 
the scintillators. Five layers of 70 //m thick Teflon ribbon are wrapped around the four 
side walls of each scintillator and two layers of single-sided 6 jum thick aluminized Mylar 
foil are wrapped around the entrance face. One layer of double-sided 12 //m thick alu- 
minized Mylar foil is wrapped around the other five faces to protect against light produced 
by Geiger discharge in the tracking region. Bicron BC600 optical glue with a nominal 
thickness of 100 /im is used to couple the light guides to the scintillator blocks. 
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Table 1 

Types and dimensions of NEMO-3 inner and outer wall scintillator blocks. The blocks 
were made at the Joint Institute for Nuclear Research (JINR) in Dubna, Russia. [[7] 



Block type 


IN 


EC 


EE 


Thickness 


98 to 110 mm 


99 mm 


99 to 123 mm 


Height 


153 mm 


200 mm 


200 mm 


Width 


138 to 154 mm 


218 mm 


218 to 230 mm 


Associated PMT [Q2] 


R6091 (3 inch) 


R6594 (5 inch) 


R6594 (5 inch) 


Total number 


680 


260 


520 



2.2. Input components of the simulation model 

We use spectral properties of all optical elements. The effective quantum efficiency 
spectraH for the PMTs are shown in Figure [2j where data are taken from [ [12] . The 
refractive index data for borosilicate glass are given by the Cauchy dispersion law n g i ass = 
1.472 + 3760/ A 2 , where A is the photon wavelength [[13]. Figure [2] also shows the refractive 
index input data for the PMMA light guides and for scintillator polystyrene [Ej, which 
we have linearly extrapolated in the low wavelength (200 — 300 nm) region. 

The measured reflection coefficients of Teflon and aluminized Mylar used in the simu- 
lations are shown in Figure [3j We assume a 50 jum air gap between the scintillator block 
and the Teflon/ Mylar wrappings whose reflectivities are modeled as 100% Lambertian 
for Teflon and 100% specular for aluminized Mylar. The scintillator blocks were treated 
under water with 1200 grit sandpaper to obtain diffusive reflection at the surfaces [[?]. To 
address this, we incorporate the surface roughness parameter "sigma_alpha" of GEANT4 
with a value of 360 degrees. 

We treat the polystyrene and the primary dopant, pTP as a single entity and define 
the primary emission spectrum to be that of the pTP alone, since the number of emitted 
photons from the polystyrene is negligible compared to the number of photons generated 
by pTP. The pTP emission spectrum is shown in Figure [H Photons emitted from pTP can 
either propagate in the polystyrene according to the bulk absorption length (BAL) of the 
polystyrene and pTP mixture or interact with a POPOP wavelength-shifting molecule. 
Stokes shifting is determined by the combined absorption length of PS/pTP/POPOP 
with effective emission governed by the POPOP emission spectrum using data shown in 
Figures [3] and [H We account for the molecular quantum yield of POPOP by having 
absorbed photons re-emitted with 93% probability [ [15] and at a wavelength equal to or 
greater than the absorbing wavelengthjf to obey energy conservation [US]. In GEANT4, 
two BALs must be specified: one for processes where the primary photon is absorbed in 
PS/pTP, and one for processes where the primary photon is absorbed by POPOP with 
the possibility of multiple wavelength shifts. 



2 Effective quantum efficiency is the product of quantum efficiency and coilection efficiency. 

3 A more precise treatment is possible for example, using Jablonsky diagrams for energy level spacing but 

would require an additional layer of complexity which we chose not to introduce at this time [ [16] . 
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Figure 1. One of twenty sectors of the NEMO-3 detector with details showing the source 
foil, scintillator blocks, and photomultipliers. EE, EC, and IN identify blocks on the 
exterior and interior walls. LI - L4 identify blocks on the petals (not modeled in this 
work). The lower figure shows a 2-D rendering of the wall blocks. 



3. Energy calibration and impact corrections 

For energy calibration and alignment of blocks, the NEMO-3 experiment uses radioac- 
tive sources that are periodically inserted into the apparatus in well-defined positions. 
Measurements are taken using sixty 207 Bi sources with a mean activity of approximately 
210 Bq. Three sources per sector are used and each twenty- four hour run yields approxi- 
mately two (three) thousand useful electron tracks for each IN (EE and EC) scintillator 
module. The bismuth decay provides conversion electrons with energies of 482 keV, 
976 keV, and 1682 keV (K-lines). The electrons lose energy due to crossing several ma- 
terials including the kapton window of the calibration tube [ [7J , the helium gas, and the 
scintillator wrapping. The mean energy losses are estimated to be 45 keV and 40 keV for 
the 482 keV and 976 keV electrons, respectively. The relation between the PMT charge 
signal and the energy deposited in the scintillator block is linear from 150 keV to 4 MeV. 
We apply a linear formula to calculate the energy of each electron event: E = a(C — P) + /3 
where C is the ADC value of the scintillator, P is the pedestal and a and are fit param- 
eters. The fit takes into account the energy loss calculated for each electron according to 
its energy and measured track length. The energy resolution, Rfwhm at 1 MeV, for each 
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Figure 2. Input data used in the NEMO-3 simulations. Left: Typical quantum efficiency 
of Hamamatsu R609I 3 inch and Hamamatsu R6594 5 inch PMTs. Right: Refractive 
indices of borosilicate glass, polystyrene, and PMMA. 
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Figure 3. Other input data used in the NEMO-3 simulations. Left: Reflection coeffi- 
cients for the Teflon tape and aluminized Mylar. Right: Absorption lengths for PMMA, 
polystyrene, and the combined polystyrene/pTP/POPOP scintillator, where the absorp- 
tion data were taken from [ [5] and [ [17] . The combination defines the scintillator absorp- 
tion length after introduction of the fluors to the polystyrene. 
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Figure 4. Primary emission spectrum for pTP and POPOP used in the simulations [[15]. 



block, is obtained from the width of the 976 keV peak assuming 
_ FWHM(E) 

iv-FWHM^J — ™ 

where A(MeV) is a constant with a range from 0.014 to 0.032 for all blocks. The mean 
value of Rfwhm at 1 MeV for EC blocks is 13.8%, for EE blocks is 13.5%, and for IN 
blocks is 16.7%, as summarized in Table |2j 

The response of each scintillator block depends upon the entrance point of the electron. 
NEMO-3 data show a weak dependence of 1 to 2% for blocks equipped with 3 inch PMTs 
and a stronger dependence of up to 10% for larger blocks equipped with 5 inch PMTs. 
This effect has a non-negligible consequence on the energy resolution. The ADC value 
and energy loss were measured for each electron in each scintillator block to correct for 
this dependency. The front faces of the scintillator blocks were divided into nine (3 x 3) 
rectangles of equal areas for modules equipped with 3 inch PMTs and twenty-five (5 x 5) 
rectangles of equal areas for modules equipped with 5 inch PMTs. ADC histograms 
were obtained for each impact region. Then the energy distribution around the 976 keV 
electron peak was refit to obtain the parameter A of the resolution function. The mean 
energy resolutions Rfwhm (E) at 1 MeV calculated using these corrections are given in 
Table H 

4. Results of simulations 

The goal of our simulations was to understand the measured energy resolution and 
the response non-uniformity as a function of impact position. We investigated the light 
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Table 2 

Simulated and measured mean energy resolutions at 1 MeV for three different types of 
NEMO-3 scintillator modules. Simulation uncertainty is taken from Table [3]for 10% input 
variations. Measurement uncertainty is calculated from the collective response variability 
between individual blocks. 



Block Type 


Simulated Rfwhm 


Measured Rfwhm 
with impact corrections 


Measured Rfwhm 
without impact corrections 


EC 


14.4 ± 1.1 


13.8 ± 0.3 


15.6 ± 0.3 


EE 


14.0 ± 1.1 


13.5 ± 0.2 


14.9 ± 0.2 


IN 


14.9 ± 1.1 


16.7 ± 0.2 


16.8 ± 0.2 



collection from electrons incident on the faces of three types of NEMO-3 scintillator blocks. 
In our model, 1 MeV electrons were generated 60 cm away from each type of scintillator 
block in vacuum, for direct comparison with the NEMO-3 geometry and the measured 
energy resolution. The angular distribution of simulated electrons matched the solid angle 
of the region of interest on the block surface. The block faces were divided in the same 
fashion as the measured data. We simulated 2,500 electrons with energy of 1 MeV in each 
region of the grid. Following the literature, we assumed the light yield for the polystyrene- 
based NEMO-3 blocks to be 8, 000 photons per MeV for electrons [[5]. This value is not 
well known for our scintillator and is a source of systematic uncertainty. The simulated 
energy resolutions in three types of blocks are compared to measured values shown in 
Tabled The energy resolution Rfwhm is calculated as Rfwhm — 2.35/ ^/Np e , where N pe 
is the number of photoelectrons registered by the PMT per each simulated electron. 

4.1. Spatial non-uniformity of the response 

We compared the simulated and measured spatial response non-uniformity by normal- 
izing the number of simulated photoelectrons and the measured PMT charge from each 
region to the average over the entire block, shown in Table [2j The EC block response 
uniformity is shown in Figure |^ a) and (b). The response reflects the symmetry of the 
block itself. For the EE block, one side is deeper than the other which results in a skewed 
response as is clearly seen in Figure |6]^a) and (b). The smaller IN block displays better 
uniformity although, one can still clearly identify the structural shape of the block in the 
results shown in Figure [TJ^a) and (b). Our simulations reproduce the measured response 
non-uniformity of each NEMO-3 block type. The number of simulated photoelectrons 
collected by each type of module is shown in Figures EJc), E(c), and He). The ratio of 
the normalized simulation value to the normalized measured value show variations no 
greater than 2% for all block types, as shown in Figures [5]^d), |6^d), and[7]^d). The ratio 
of the minimum to maximum response in collected photoelectrons is 86% and 92% for 
EC and EE simulations, respectively, and 88% for measurement in both block types. The 
ratio of the minimum to maximum response is 97% for IN block simulations and 98% for 
measurement. The statistical uncertainty in the mean number of photoelectrons collected 
for each bin is 0.04%. 
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(c) Profile of Detected Photons (d) Ratio = Simulation/Measurement 




Figure 5. Results from 1 MeV electrons incident on the EC block. Simulated (a) and 
measured (b) response, normalized to the mean response, (c) The mean number of photo- 
electrons collected in each sub-region, (d) The ratio of simulation to measurement. 



4.2. Characteristics of simulated photons 

Figure M shows the incident angle of detected photons measured from the normal to the 
photocathode surface for the EC block. FigureElalso shows the initial and final wavelength 
profiles of the detected photons for the EC block. The distribution resembles that of 
the POPOP emission spectrum. Shorter wavelength photons are absorbed in the bulk 
scintillator material and on the surface of the wrappings resulting in a suppression below 
400 nm. Longer wavelength photons are supressed due to the low PMT quantum efficiency 
at longer wavelengths. Figure [H] also shows the detection probability as a function of the 
number of wavelength shifts per photon, the number of reflections from the wrapping 
surfaces, and the number of total internal reflections. Each individual photon, for example, 
can undergo several wavelength shifts and several reflections (specular, diffuse, and total 
internal) before reaching the PMT. The distribution of wavelength shifting events suggests 
that a large fraction of the photons from pTP have been wavelength shifted. However, 
the probability for zero wavelength shifts is finite because pTP emission can take place 
above 400 nm, out of the range of large bulk attenuation. 

4.3. Dependence of simulations on input parameters 

For the EC block, we have varied the simulation input parameters to evaluate the 
dependence of the energy resolution on these changes. The PMT quantum efficiency, the 
absorption length of the scintillator, the reflectivity of Teflon and Mylar wrappings, and 
the light yield were decreased and increased by 5% and 10%. These fractional changes 
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(a) Simulation 



(b) Measurement 




(c) Profile of Detected Photons 



(d) Ratio = Simulation/Measurement 
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Figure 6. Results from 1 MeV electrons incident on the EE block. Simulated (a) and 
measured (b) response, normalized to the mean response, (c) The mean number of photo- 
electrons collected in each sub- region, (d) The ratio of simulation to measurement. 



reasonably reflect physical variations in production or manufacturing of most components. 
We then compared the result of each change to the central value of the energy resolution 
of 14.4%. Table [3] summarizes our studies. 

4.4. Dependence of energy resolution on spectral properties 

We evaluated the importance of using spectral properties of materials in our simulations 
of energy resolution of the EC block. We began with a simplified simulation of optical 
materials with wavelength-independent values and investigated the effect of systemati- 
cally introducing the wavelength dependence into the simulation. Initially, we fixed the 
value for the quantum efficiency at 21%, the absorption length at 4.5 m, the reflectivity 
for the Teflon and Mylar at 93%, and the refractive indices at 1.5. The wavelength- 
dependent quantum efficiency of the PMT was introduced first. We then introduce a 
wavelength-dependent absorption length for the scintillator with no wavelength shifting 
and subsequently introduce wavelength-shifting dependence and molecular quantum effi- 
ciency Finally, we introduced wavelength-dependence for the Teflon and Mylar reflection 
coefficients, and the scintillator and borosilicate glass refractive indices. The values here 
were chosen as best-guess estimates at a peak wavelength of 420 nm that one may assume 
for a monochromatic simulation. The results for each step are shown in Table HI The 
complete model reproduces the measured energy resolution within systematic uncertainty. 
It has been suggested that there is a variation in the PMT quantum efficiency as a func- 
tion of the incident angle of the photon [US]- If we include this effect, assuming it is also 
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(c) Profile of Detected Photons (d) Ratio = Simulation/Measurement 




Figure 7. Results from 1 MeV electrons incident on the IN block. Simulated (a) and 
measured (b) response, normalized to the mean response, (c) The mean number of photo- 
electrons collected in each sub-region, (d) The ratio of simulation to measurement. 



relevant for NEMO-3 PMTs, the energy resolution for the complete model improves by 
about 0.3% while the spatial response distribution is left unchanged. 



5. Modeling of SuperNEMO scintillator blocks 

We have used our simulations to facilitate a scintillator choice for the SuperNEMO 
experiment [[21]. SuperNEMO further plans to exploit the NEMO-3 technique of tracking 
and calorimetry. The new modular detector would incorporate about 100 kg of 82 Se, 
150 Nd, or 48 Ca, to reach neutrinoless double beta decay half-life sensitivity of about 1.5 x 
10 26 years. This goal requires that the new experiment significantly improves its energy 
resolution with respect to NEMO-3. We have been conducting R&D and the developed 
simulation code is an important aid in our studies of the choice of materials, block shape 
and size, wrapping, and light collection. 

A baseline provisional design for the SuperNEMO calorimeter calls for a hexagonal 
block with a proposed circumscribed radius of 22.5 cm made out of polyvinyltoluene 
(PVT) scintillator (e.g., Eljen EJ-200 [ 18J) coupled to a super-bialkali 8 inch hemispher- 
ical PMT (e.g., Hamamatsu R5912-MOD). Figure M shows a drawing of the proposed 
SuperNEMO scintillator module including the scintillator, PMT, and mounting brackets. 
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Figure 8. Left: Incident angle of detected photons relative to the photocathode normal 
for the EC block. Center: Initial and final wavelengths of detected photons simulated for 
the EC block. The initial wavelength distribution is primarily that of the pTP emission. 
The final wavelength distribution is reminiscent of the POPOP emission spectrum with 
degradation due to absorption at short wavelengths and poor PMT quantum efficiency 
at long wavelengths. Right: The detection probability per photon as a function of the 
number of wavelength shifts, number of reflections off the wrapping surfaces, and number 
of total internal reflections. 



The light yield of the PVT scintillator is taken to be 10,000 photons per MeV [118^1 The 
quantum efficiency of the PMT was taken to be approximately 33% at 420 nm [ 12^|. 
We propose a high reflectivity aluminized Mylar from ReflecTech [ [19] around the sides 
and entrance face of the scintillator block and a Teflon wrapping on the top face near the 
PMT. Similar to NEMO-3, we assume an air gap of 50 /zm between all block faces and the 
Mylar wrapping. Our simulations show that with this configuration, assuming NEMO-3 
scintillator absorption and emission, a resolution of 7.5 ± 0.5% (FWHM) at 1 MeV [ 120] 
is expected. Using the symmetry of the block, each sixth of the hexagonal face is divided 
into 16 regions. As shown in Figure [101 the energy resolution is fairly uniform across the 
face of the block. The mean energy resolution is 7.19% (FWHM) at 1 MeV and the min- 
imum and maximum is 7.14% and 7.24%, respectively. Recently conducted preliminary 
measurements confirm our predictions [l2"Tj. 



4 This is compared to a polystyrene scintillator with nominal light yield of 8,000 photons per MeV 
5 This is compared to a 5 inch NEMO-3 PMT with a quantum efficiency of 25% at 420 nm 
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Figure 9. A drawing of a provisional SuperNEMO scintillator module showing the scin- 
tillator, PMT, and mounting brackets. 




Figure 10. Spatial dependence of energy resolution (FWHM at 1 MeV) for one sixth of a 
SuperNEMO scintillator block. 
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Table 3 

Changes of the EC block energy resolution at 1 MeV due to variations of input parameters. 
We increased and decreased the parameter values by 5% and 10% and quote the difference 
in the energy resolution between the central value of 14.4% and an average of the change. 
Total change is calculated by adding the individual changes in quadrature. 



Optical Parameter 


5% 


10% 


Quantum Efficiency 


0.34% 


0.75% 


Absorption Length 


0.12% 


0.31% 


Teflon Reflectivity 


0.92% 


2.50% 


Mylar Reflectivity 


0.54% 


0.98% 


Light yield 


0.30% 


0.75% 


Total change 


1.17% 


2.90% 



Table 4 

Simulated FWHM of energy resolution at 1 MeV of EC scintillator blocks in NEMO-3 for 
increasingly comprehensive model parameters. The initial values are best-guess estimates 
at a peak wavelength of 420 nm that one may assume for a simplified monochromatic 
simulation. Wavelength-dependence is introduced systematically to show the effect. 





Rfwhm 


All quantities fixed 




quantum efficiency : 21% 




absorption length : 4.5 m 




Teflon/Mylar reflectivity : 93% 




refractive indices : 1.5 


11.2 


After introduction of A dependence to: 




Quantum efficiency (at constant abs. length) 


11.6 


Absorption length 


13.8 


Absorption length with Stokes shifting 


14.2 


Mylar reflection coefficient 


14.4 


Teflon reflection coefficient 


14.4 


All refractive indices 


14.4 


Complete model 


14.4 


Measured 


13.8 
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6. Conclusions 

We have constructed a GEANT4-based spectral model of NEMO-3 and SuperNEMO 
calorimeter blocks and compared the results of simulations with measurements of single 
electrons from 207 Bi sources. The measured energy resolution and spatial dependence for 
the blocks were demonstrated to be in good agreement with simulations. The EC blocks 
show a doubly-symmetric distribution of response about the center of the block while 
the EE blocks show a symmetric distribution about the z-axis, both in agreement with 
observations. For the SuperNEMO block, the main elements improving the resolution 
over the NEMO-3 blocks are the 25% increase in scintillator light yield, the 35% increase 
in the PMT quantum efficiency, the larger PMT size, and the incorporation of the high 
reflectivity aluminized Mylar. Additionally, by directly coupling the scintillator and PMT, 
we render light guides unnecessary which further improves the transparency and response 
uniformity. 

As expected, including the spectral properties of all materials and incorporating wavelength- 
shifting absorption and emission with the effects of fluorescent quantum yield in the scintil- 
lator improves the agreement of modeling with measurements. Our approach is necessary 
for a detailed understanding of high resolution plastic calorimeters. The importance of 
these spectral simulations increases with distances travelled by photons so such models 
are necessary to study large-scale scintillator detectors. 
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